Variable range hopping in thin film with large dielectric constant 



OO 
O 

o 

U 

Oh! 
< 



c 

I 

CZ3 



C 

O 

> 

m 
m 
cn 

o 

oo 
o 



X 



B. I. Shklovskii 

Theoretical Physics Institute, University of Minnesota, 

(Dated: April 7, 2008) 



Minneapolis, Minnesota 55455 



In a film with large dielectric constant re the electric field of an electron spreads inside the film 
before exiting the film at large distances of order of red (d is the film width). This leads to the 
logarithmic Coulomb repulsion between electrons and modifies the shape of the Coulomb gap in 
the density of localized states in a doped film. As a result the variable range hopping conductivity 
in such a film has a peculiar temperature dependence, where the domain of the lncr(T) oc (To/T) p 
dependence, with the index p ~ 0.7, is sandwiched between the two domains with p = 1/2. 



Variable range hopping (VRH) is the generic mecha- 
nism of the low temperature transport in systems with 
localized electron states. When electrons repel each other 
via the Coulomb potential energy 



V(r) = -, 
Kr 



(1) 



where re is the dielectric constant of the solid the density 
of localized states g(E) has the soft Coulomb gap. As a 
result the VRH conductivity a obeys the Efros-Shklovskii 
(ES) law 0, [| 



er(T) = ctq exp 



1/2' 



(2) 



both in two-dimensional (2D) and three-dimensional 
(3D) cases. Here the characteristic temperature 



T n = 



Ce 2 



HO 



(3) 



(we use the energy units for the temperature T) , a is the 
localization length of electrons and C is the numerical 
coefficient 0j. 

This paper deals with the situation when the Coulomb 
interaction between electrons has a more complicated 
form. We consider a thin film with the thickness d> a 
and the large dielectric constant re ^> 1. The film is 
surrounded by the media with much smaller dielectric 
constant, for example just by the air with n ext = 1. 

The energy of the Coulomb repulsion of two electrons 
in the film was calculated exactly in Ref. Here we 
present only asymptotic results with their physics inter- 
pretation. Let us assume that the film is defined by the 
surfaces z = ±d/2 and one electron is at z = a; = y = 0. 
Then at distances r <C d its electric field (induction) 
spreads isotropically. At larger distances the electric field 
lines are forced by the large re to stay inside the film, so 
that the field spreads along the radius p = \J x 1 + y 2 of 
the cylindrical coordinate system with the same z-axis. 
At p ~ Hid electric field lines exit from the film and even- 
tually spread uniformly over the whole 4-7T body angle 
again. 

Let us discuss the potential energy of repulsion of two 
electrons in the film. At the distance p 3> red two elec- 



trons interact via the "external" Coulomb interaction 

V(p) = -. (4) 
P 

On the other hand, in intermediate range of distances 
d <C p <C red 



V(p) = 2,,\u(jj 



(5) 



where = e 2 / 'red. Finally at even smaller distances p <C 
d we arrive at the the "internal" Coulomb interaction 



V(r) = — + 2e d lnre, 



(6) 



which differs from Eq. (Q]) by the logarithmically large 
energy accumulated when the second electron is moved 
from infinity to p = d. The resulting potential energy 
V(p) is plotted as a function of p in Fig. Q] in all three 
ranges. 



V(p)f 



(6) 



(5) 



(4) 



d 



Kd 



FIG. I: A schematic plot of the potential energy of repulsion 
of two charges located inside the film with a large dielectric 
constant as a function of the distance p between electrons. 
Equation numbers describing different segments of this plot 
are shown next to each segment. 

Let us assume that the film is uniformly doped by a 
concentration of donors No and compensated by smaller 
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concentration of acceptors so that low temperature trans- 
port is due to VRH on donors. How does the two- 
dimensional electrostatics affect ES law? 

This question was first addressed in Ref. [J] for an- 
other object, a weakly disordered two-dimensional gas, 
for example, in a silicon inversion layer. In this case, ac- 
cording to Ref. [f| localized states have an exponentially 
large localization length. Large localized states play the 
dual role in this theory. First, VRH conductivity is re- 
lated to long distance hopping between those of them 
which are close to the Fermi level. Second, all other large 
localized states according to Ref. [|| contribute to the 
large effective dielectric constant which keeps the elec- 
tric lines of a charge in the plane of the two-dimensional 
gas. In order to calculate VRH conductivity the authors 
of Ref. Tj used an intuitive shortcut avoiding discussion 
of the Coulomb gap. The authors found two ranges of the 
temperature dependence of the VRH conductivity a(T). 
At very low temperatures where the characteristic length 
of the hop r% is much larger than nd they arrived at ES 
law with k = 1 and and at higher temperature range 
where d <C <C nd they obtained activated conductiv- 
ity. 

In this paper, we calculate the VRH conductivity in 
the framework of our simpler model, where large k is of 
the lattice origin or is a result of close three-dimensional 
metal-insulator transition. In this way we avoid the con- 
troversial subject of the 2D metal-insulator transition Q. 
We follow the "orthodox" ES logic P, 0] deriving first 
the Coulomb gap and then the conductivity. Let us first 
formulate our results moving from high to lower temper- 
atures. 

At high temperatures T 3> e^a/d the VRH hopping 
conductivity is determined by the "internal" ES law 



the power 



er(T) = (To exp 



1/2 



£d In K 
T 



(7) 



which at high enough temperatures coincides with 
Eq. ©. 

It is followed by the range of intermediate tempera- 
tures 



e -^« T « tJ T 
nd d 

with "activated" VRH conductivity 

Ti(T) 



cr(T) = do exp 
where the "activation energy" 
T 1 (T)=e d In 



T 



T nd 
td a 



(8) 



(9) 



(10) 



of the intermediate regime is weakly temperature depen- 
dent. If one approximates Eq. (J9j) as 



cr(T) = (Tq exp 



1 p 
T 



p=l 



hi{TK,d/eda) 



(12) 



Close to T = e^a/d we get p = 1 — 1/ln/c. For example, 
at k — 40 one gets p — 0.7. 

At even smaller temperatures T <C erfo/fcci we arrive 
at the "external" ES law: 



<j(T) = (To exp 



rpext 
±0 

T 



1/2' 



(13) 



(11) 



where T§ xt = Ce 2 /a. Thus, the "activated" VRH con- 
ductivity Eq. © is sandwiched between the two different 
ES regimes of er(T), the "internal" one, Eq. 0, at the 
high temperature side, and the "external" one, Eq. (flUj) , 
on the low temperature side. The two low tempera- 
ture regimes Eq. (|9|) and Eq. () 1 3(1 are in agreement with 
Ref. [4J. The new high temperature "internal" ES regime 
exists only if d 3> a. 

The experimental literature on the VRH conductivity 
in thin films is controversial (see Refs. @, Q and ref- 
erences therein.) How large is the film dielectric con- 
stant and how important is contribution of large localized 
states [4( in most of cases is not clear. On the insulating 
side of the superconductor-insulator transition in ultra- 
thin quench-condensed Ag, Bi, Pb and Pd films [7j VRH 
data agree with Eq. (TTTj) with p = 2/3. On the other 
hand, for relatively thick TiN films [|| the crossover from 
p ~ 1/2 to p ~ 1 is observed with the decreasing tem- 
perature. Finding an explanation for this crossover is 
challenging [8j because one would normally expect that 
the ES law emerges at the low temperature limit. Our 
theory shows that in relatively thick film with the de- 
creasing temperature one can see the crossover from the 
"internal" ES law to the activated transport. This may 
explain results of Ref. @ . 

In order to simulate a film with a large dielectric con- 
stant one can make a two-dimensional array of isolated 
metallic islands overhanging each other Although 
these arrays were originally designed as arrays of Joseph- 
son junctions they perfectly simulate a large dielectric 
constant in the normal state. Indeed, such array keeps 
electric lines in its plane if the capacitance between two 
islands is larger than the capacitance of each island to 
the ground. As a result such a normal array in presence 
of some disorder should show activated VRH discussed 
above. 

Before switching to the derivation of our results we 
would like to dwell on the related theoretical paper [ll[ 
which deals with the VRH transport of point like vor- 
texes responsible for the low temperature resistance of a 
superconductor film in the external magnetic film. Two 
vortexes interact via the logarithmic potential at small 
distances, while again at large distances their interaction 
follows the "external" Coulomb potential [l|, so that 
one could expect to see the two low temperature ranges 
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discussed above, the "activated" regime and the "exter- 
nal" ES law. However, Ref. [ll[ argues that ES approach 
is not valid in this case because "logarithmic interac- 
tion grows without bound with particle separation" and, 
therefore, "single-particle energies can not be defined" so 
that " multi- vortex hopping dominates the above single- 
particle effects". The multi-particle estimate (TlJ leads 
to Eq. (JTTJ) with 2/3 < p < 4/5. This is close to what we 
got for p at k — 40. 

In contrary to the above statements of Ref. Fig. [1] 
clearly shows that the repulsion energy vanishes at infin- 
ity similarly to the standard Coulomb potential Eq. ([J|. 
Thus, there is no problem to introduce a single particle 
energy for a system of localized electrons interacting with 
the pairwise potential V(p). We can proceed with the ES 
argument and study the new shape of the Coulomb gap 
in the density of states (DOS) of single particle excita- 
tions and eventually the VRH conductivity. We return 
to the discussion of the role of multi-particle processes in 
the end of this paper. 

Below we calculate the zero temperature DOS g(e) fol- 
lowing Refs. [l|, In the ground state of the system 
we define the single electron energy q of an occupied 
donor i as the energy necessary to extract an electron 
from this state to the state with the energy right at the 
Fermi level at infinity. The single electron energy ej of 
an empty in the ground state donor j is defined as the 
energy necessary to bring to it an electron from a state at 
infinity with the energy right at the Fermi level ep. By 
the definition of the ground state ej > ep > Another 
stronger stability condition can be formulated for each 
pair of occupied and empty states as follows 



> 0. 



(14) 



Here V(rij) is the repulsion energy of two electrons on 
sites i and j and —Vfcj) is the Coulomb energy of at- 
traction of the electron which moved to the site j by the 
hole it has left at the site i. In other words, this term 
describes the exciton effect. Eq. (|14p requires that any 
two states close in energy to ep should be far enough in 
the space. This limits the density of states (DOS) close 
to the Fermi level. For the Coulomb potential the result 
for the DOS is known and we do not repeat derivations 
from [U [2l. [l3T] . but list the results. 

The "external" Coulomb potential Eq. limits two- 
dimensional density of states at the level 



, , 2 e 



7r e* 



(15) 



For the "internal" Coulomb potential Eq. 

3 K 3 e 2 d 
■9( e ) = ! 



we get 



(16) 



where the factor d converts the three-dimensional DOS 
to the two-dimensional one. Apparently the asymptote 
Eq. (fT5"l) is valid at large distances, i. e. at small energies 
e and the asymptote Eq. (fT5|) describes small lengths or 



large energies. Let us now derive the DOS for the in- 
termediate range of distances and energies. To this end 
we have to calculate the number of states in the band of 
the width 2e around ep. Using Eq. (fT4"|) and Eq. ([5]) we 
get ed \n{nd/ p) < e. This means that in this energy band 
there is no more than one state in the disc with radius 



p(e) = ndexp 



e 



(17) 



This leads to the following estimate of the DOS for td <C 
e -C e d ln/t: 



9(e) 



1 



1 



p(e) 2 e {dnf 



■ exp 



(18) 



DOS Eq. HH1) matches DOS Eq. (JI5J) at e = e d and 
Eq. (fT6|) at e = e^lnK (see Fig. [2])). At large energies 
the parabolic range of the Coulomb gap, Eq. (fTB")) . is lim- 
ited by the total width of the impurity band e max = 
e 2 Np 3 /n + e^lnK. At £ > (-max the DOS decreases 
(see the dashed line in Fig. [5]) similarly to the three- 
dimensional DOS of the classical impurity band (see 
Ch.14 of Ref. 0). 




FIG. 2: Schematic plot of the DOS of localized electrons lim- 
ited by Coulomb interaction of electrons as a function of en- 
ergy. Equation numbers describing different segments of this 
plot are shown next to each segment. Dashed line shows be- 
ginning of the DOS decline at large energies 

Now we can apply the calculated above DOS in order 
to estimate the exponential term of the VRH conductiv- 
ity of the film at different temperatures. This estimate 
closely follows the original Mott's approach [14} . We de- 
fine an energy band around the Fermi level and estimate 
contribution of this band to the VRH conductivity. Then 
we can optimize result with respect of the band width. 
This calculations are quite straightforward and we will 
not go through them here. Results are already formu- 
lated in the beginning of our paper. 

Let us return to the multi-electron effects on ES law. 
For the standard Coulomb potential Eq. |T]) they were 
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studied in Refs. @, [H, [IH, [l6| • It was shown that they 
may change only the coefficient C in Eq. @. Here 
we only briefly remind what was done. In 3D case a 
small energy single-electron excitation strongly interacts 
with dipole moments of surrounding compact electron- 
hole excitations forming together with them a composite 
charged multi-electron excitation, the electronic polaron. 
Polarons being charged particles obey stability criterion 
Eq. ([M]) , have the Coulomb gap and lead to ES law. The 
only difference is that coefficient C can be somewhat 
larger because every hop of a single-electron excitation 
in 3D is accompanied by the tunnelling depolarization 
of many polaron pairs. A simple estimate showed that 
the total length of these small hops is of the order of the 
length of the main long ES VRH hop. This gives no more 
than the factor 2 in the expression for C . In 2D number 
of dipole excitations in the polaron atmosphere is of order 
one and polaron effects provide only small corrections to 
C. 

We estimated similar effects in the framework of the 
complex potential of this paper (Eqs. ([JJ, ©). It 
turns out that in the low temperature range polaron ef- 
fects provide only small corrections to C in Eq. (|13p . On 
the other hand, for the two higher temperature ranges 
(Eqs. (JOj) and ((2]) the situation is similar to the 3D case 
studied in Refs. @, [rJ [l6| . Namely, at these tempera- 
tures multi-electron effects may add a numerical coeffi- 
cient in the exponential. For Eq. ([9]) this means that T\ 
may become twice larger. Thus, multi-electron effects do 
not change the power p in Eq . (jlip. This conclusion is 
in disagreement with Ref. [lll | which assumes that even 
in 2D the total length of small hops is much larger than 
the ES hop and, therefore, overestimates importance of 
the multi-electron effects. 

Above we discussed only the ohmic transport in a 
weak electric field. If the electric field is so strong that 
eEa/2 ^> T one can replace T by the effective VRH 
temperature eEa/2 in Eqs. ©, ©, ([11]) and |T31) to ob- 
tain non-ohmic current-voltage characteristics [17| (see 
also [3, For the intermediate activated regime we 

then arrive at the current-voltage characteristics J oc 
Joex-p(—Eo/E) p , where p is close to unity. This result 



is in agreement with the earlier theory [l9[ of the VRH 
transport of pinned vortexes in superconductors under 
the influence of a strong current. 

Until now we dealt with a film. Now let us briefly dis- 
cuss the similar physics in a long cylindrical nano-wire 
or nano-rod with radius d made from a semiconductor 
with a large dielectric constant k>1, for simplicity, in 
the air environment. Here again the electric field of an 
electron at distances r <C d spreads isotropically, then 
stays inside nano-rod for a distance x < £ — dn 1 / 2 along 
the cylinder axis, and then leaks from the cylinder and 
eventually spreads isotropically in the air at large enough 
r. Thus, the potential energy of repulsion of two elec- 
trons again changes from the very short range "inter- 
nal" Coulomb interaction Eq. (j4|) to the very long range 
"external" Coulomb interaction Eq. ((T|) with the large 
intermediate range of x, where interaction has the one- 
dimensional character. The interaction in this range was 
studied [20 . |21| for an ion channel in a lipid membrane, 
where the cylindrical pore with radius d is filled by wa- 
ter with k = 81 and is surrounded by lipids with n = 2. 
Translated to our problem the potential energy of two 
electrons located at the nano-rod axis at the intermedi- 
ate distance x from each other is well approximated by 

V(x)=eEa{[exp(-x/Z)-l], (19) 

where E = 2e/nd 2 . (see Sec. VIII of Ref. (HJ). Using 
this potential for repulsion of two electrons in the nano- 
rod and following the Coulomb gap based derivation sim- 
ilar to one used above for a film (or the shortcut approach 
of Ref. L 4]) one can calculate the temperature dependence 
of the VRH conductivity of the nano-rod in the interme- 
diate temperature range. The result is the strict activa- 
tion regime with the temperature independent activation 
energy T a — eE^. It is, of course, sandwiched between 
the two ES laws, the "internal" one on the high temper- 
ature side and the "external" one Eq. ([i~3")) on the low 
temperature side. 

I am grateful to T. Baturina, A. M. Goldman, A. 
Kamenev, D. E. Khmelnitskii and A. I. Larkin for useful 
discussions. 
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